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Q ' Abstract. Particle induced X-ray emission (PIXE) is a physical effect that is not yet 

O , adequately modelled in Geant4. The current status as in Geant4 9.2 release is reviewed and 

' new developments are described. The capabilities of the software prototype are illustrated in 

O ' application to the shielding of the X-ray detectors of the eROSITA telescope on the upcoming 

' ^ , Spectrum-X-Gamma space mission. 

>>■ 

^ : 

y—^ , 1. Introduction 

' The wide use of PIXE as an experimental technique has motivated several dedicated software 

■ systems; nevertheless, limited functionality for PIXE simulation is available in general-purpose 
I Monte Carlo codes. Geant4 [H [2] addresses X-ray emission induced both by electrons and heavy 

fN| ' particles like protons and a particles. 

. The simulation of PIXE involves the energy loss and scattering of the incident particle, atomic 

I shell ionization cross sections, and atomic transition probabilities and energies. Intrinsically, 

■ PIXE is a discrete process: X-ray emission occurs as the result of a vacancy in shell occupancy, 
. . i in competition with Auger electron emission and Coster-Kronig transitions. Nevertheless, this 

' discrete process is intertwined with the ionization process, which determines the production of 

^ . the vacancy; this process is treated in general-purpose Monte Carlo codes with mixed condensed 

^ I and discrete transport schemes, since excessive computational resources would be necessary to 

■ ■ ' handle the infra-red divergent cross section with a discrete scheme. 

While this combined condensed and discrete particle transport scheme is appropriate to 
many simulation applications, it suffers from drawbacks with respect to the generation of PIXE. 
Atomic relaxation occurs only in connection with the discrete part of the transport scheme, 
where the production of a 5-iay can be associated with the creation of a vacancy; therefore 
the fluorescence yield depends on the threshold for the secondary electron production. Another 
drawback is that the cross section for the production of a 6-ray is calculated from a model for 
energy loss that is independent of the shell where the ionization occurs. 

This paper analyzes tools pertinent to PIXE simulation available in Geant4 at the time of the 
CHEP 2009 conference and describes a set of developments to extend the simulation capabilities 
in this domain. A more extensive account can be found in |3|. 
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Figure 1. A comparison of K shell ionization cross sections, implemented in 
G4hShellCrossSectionDouble [101 HI] (black circles), with Paul and Bolik [13j empirical reference 
for a particles (open symbols): the curves are respectively for Si (squares), Cu (triangles), Cd 
(crosses) and Au (white circles). The line represents ECPSSR K shell cross sections for a 
particles from the prototype developments of this paper. 

2. Current Geant4 tools for PIXE simulation 

Software models for the simulation of PIXE induced by protons and a particles are currently 
available in the Geant4 low energy electromagnetic package [H [S] . 

2.1. First set of models 

The first development cycle [6] involved a software design for interchangeable ionization cross 
sections models. Based on the cross section calculations, the ionisation process creates a vacancy, 
and uses Geant4 Atomic relaxation [7J for producing fiuorescence photons and Auger electrons. 
Three ionization cross section models were implemented. 

The G4hSheUCrossSection class implements cross sections based on Gryzinski's model [8^, '9]. 
It appears to be affected by a severe problem of inconsistency with experimental data. 

K shell ionization cross sections based on fits to tabulated data are implemented |10l [TT] 
in G4hShellCrossSectionExp and G^hShellCrossSectionDoupleExp. These models have been 
documented by their implementers [lOt [TT] as cross sections for K shell ionization by protons 
based on Paul and Sacher compilation [12]. Nevertheless, they appear to describe the cross 
sections for K shell ionization by a particles based on by Paul and Bolik compilation [Ij], as it 
is evident in Fig. [TJ 

The algorithm implemented to identify the shell where a vacancy is created introduces 
a dependency on the user-defined electron production threshold. This is prone to generate 
inconsistent behavior, as it is illustrated in Fig. |21 A detailed discussion of this issue can be 
found in [3]. 
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Figure 2. The filled symbols are Si ionization cross sections from an ECPSSR calculation: 
K shell (black circles) and L shell summed over all sub-shells (black squares). The total cross 
section for the discrete part of ionization as computed in Geant4 is depicted with open symbols 
for different 5-iay production thresholds: 1 keV (white circles), 10 keV (white squares) and 
100 keV (white triangles). For some 6-Tay production thresholds the total cross section results 
smaller than the cross section for ionising any shell. 



2.2. Models in Geant4 9.2 

A set of ionization cross section implementations [T^ for FIXE simulation were released in 
Geant4 9.2. The cross sections for K shell ionization are based on interpolations of empirical 
data compilations (Faul and Sacher [12] for protons, Faul and Bolik [13] for a particles) and on 
the ECFSSR theory [15j; they are respectively implemented in the Geant4 G4PaulKCrossSection 
and G4ecpssrCrossSection classes. The cross sections for L shell ionization by protons are based 
on the semi-empirical model of |16j . 

The plots in [14j compare Geant4 9.2 cross section implementations for K shell ionization by 
protons and proton ionization experimental data, with the cross section models for a particles 
of the first development cycle: this comparison is not pertinent to demonstrate the capabilities 
of the Geant4 9.2 implementations with respect to the previous ones. 

The software described in pH] and released in Geant4 9.2 is affected by several drawbacks. 

G^PaulKCrossSection is slower in execution by a factor 40 to 60 approximately than the 
previously existing classes implementing the same functionality, and two orders of magnitude 
slower than the equivalent implementation described in this paper. Other drawbacks are a 
memory leak and the attempt to open data files not included in the Geant4 9.2 release. 

The G4ecpssrGrossSection class computes K shell ionization cross section for protons and 
a particles nominally based on the ECFSSR formulation |15j . However, it does not appear to 
reproduce ECFSSR cross sections correctly: for instance. Fig. 1 to 3 in |14j exhibit an evident 
divergence at higher energies of these implementations with respect to the empirical data by 
Faul and Sacher. Figure [3] shows that this discrepancy hints to a flawed software implementation 
of the ECFSSR theory in G4ecpssrGrossSection, rather than a deficiency of the ECFSSR theory 
as suggested by the authors of [Tl]. This class results in slower execution (by a factor 6-8) than 
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Figure 3. Cross sections for copper K shell ionization by protons: reference ECPSSR 
tabulations reported in [12] (stars), ECPSSR as calculated by the G^ecpssrCrossSection class 
released in Geant4 9.2 [T^ (white crosses), ECPSSR model as calculated by the software 
described in section [3] based on ISICS tabulations (solid line), and empirical model of Paul and 
Sacher [T2] (black circles). G^ecpssrCrossSection does not appear to implement the ECPSSR 
theoretical model correctly. 

the ECPSSR models described in this paper. 

The semi-empirical model by Orlic et al. [TB] implemented in G^OrlicLCrossSection has 
limited capabilities: it computes cross sections for Li^2,3 sub-shells only for elements with atomic 
number greater than 40 and is valid only for proton energies up to approximately 10 MeV. 
Moreover, as demonstrated in section HT2l it is less accurate than other available models. 

The cross section models released in Geant4 9.2 exhibit severe software design limitations, 
that prevent client code from exploiting either dynamic or static polymorphism. Due to their 
design features, these models are not practically usable for PIXE simulation applications based 
on Geant4; their implementation features - correctness, computational performance and physical 
accuracy - motivate concerns about their use in experimental studies. 

3. Prototype developments for PIXE simulation vi^ith Geant4 

The new developments for PIXE simulation involved a re-design of the software, extended physics 
capabilities, the validation against experimental data and the application to an experimental 
use case. 

The domain decomposition at the basis of PIXE simulation with Geant4 identified three main 
entities with associated responsibilities: the hadron ionization process, the creation of a vacancy 
in the shell occupancy resulting from ionisation, the deexcitation of the ionised atom with the 
associated generation of X-rays. The simulation of PIXE is the result of the collaboration of 
these entities. A class diagram in the Unified Modelling Language (UML) illustrates the main 
features of the software design in Fig. 21 
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Figure 4. The Unified Modelling Language class diagram of the developments for PIXE 
simulation described in this paper, illustrating the main features of the software. 



3.1. lonisation cross section models 

A wide choice of cross section models for K, L and M shell ionization is provided in the prototype 
software for protons and a particles. The availability of ionization cross section calculations and 
experimental data for outer shells is very limited in literature. Theoretical cross section models 
include Plane Wave Born Approximation (PWBA) and variants of the ECPSSR model: the 
original ECPSSR formulation [E], ECPSSR with United Atom correction (ECPSSR-UA) [IT], 
ECPSSR with corrections for the Dirac-Hartree-Slater nature of the K shell [18] (ECPSSR-HS), 
as well as calculations based on recent improvements to K shell cross section specific to high 
energy [H] (ECPSSR-HE). 

The cross sections have been tabulated and assembled in a data library; the values at a given 
energy are calculated by interpolation. The tabulations corresponding to theoretical calculations 
span the energy range between 10 keV and 10 GeV; empirical models are tabulated consistently 
with their energy range of validity. The adopted data-driven approach optimizes performance 
speed and offers flexibility for chosing a cross section model. 

ECPSSR tabulations have been produced using the ISICS software [201 [21], 2006 version and 
an extended version |22j including recent high energy developments. Tabulations of ECPSSR 
calculations as reported in [12] are also provided. 

Empirical cross section models for K shell ionization include the tabulations for protons 
documented in [12] and a more recent one [23]. An empirical cross section model for K shell 
ionization by a particles is based on the tabulations in |13] . Empirical models for L shell 
ionization by protons have been developed by Miyagawa et al. [21], Sow et al. [25] and Orlic et 

ai. m- 

The ISICS software allows the calculation of cross sections for heavier ions as well; therefore, 
the current PIXE simulation capabilities can be easily extended in future development cycles. 



Table 1. Percentage of test case with compatibility at confidence level CL between simulation 
models and experimental data of K shell ionization by protons 
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3.2. Generation of a vacancy 

The determination of which atomic (sub-)shell is ionised is related to its ionisation cross section 
with respect to the total cross section for ionising the target atom. However, as previously 
discussed, the condensed-random- walk scheme raises an issue as to estimating the total ionisation 
cross section at a given energy of the incident particle. 

A different algorithm has been adopted with respect to the one implemented in the first 
development cycle: the vacancy in the shell occupancy is determined based on the total cross 
section calculated by summing all the individual shell ionisation cross sections. This algorithm 
provides a correct distribution of the produced vacancies as long as ionisation cross sections can 
be calculated for all the atomic shells involved in the atomic structure of the target element. 
Since cross section models are currently available for K, L and M shells only, at the present 
status of the software this algorithm overestimates PIXE for elements whose atomic structure 
involves outer shells, because of the implicit underestimation of the total ionization cross section. 
This approach, however, provides better control on the simulation results than the algorithm 
implemented in the first development cycle. 

The production of secondary particles by the atomic relaxation of an ionized atom is delegated 
to the Atomic Relaxation component. 

4. Softvi^are validation 

The availability of a wide variety of cross section models for the first time in the same 
computational environment allowed a detailed comparative assessment of their features against 
experimental data. Due to the limited page allocation in the conference proceedings and 
copyright constraints with the publisher of [3], only a brief summary is reported here; further 
details can be found in [3], including the full list of elements subjects to test against experimental 
data and a large collection of plots demonstrating the comparison between simulation and 
experimental data. 

The comparison of cross sections as a function of energy was performed for each element by 
means of the test. Contingency tables were built on the basis of the outcome of the t^st 
to determine the equivalent, or different behavior of model categories. The input to contingency 
tables derived from the results of the test: they were classified respectively as "pass" or "fail" 
according to whether the corresponding p- value was consistent with a 95% confidence level. The 
contingency tables were analyzed with Fisher exact test [26j . 

4.1. K shell ionization cross sections 

The reference experimental data were extracted from [12| . 

The fraction of test cases for which the test fails to reject the null hypothesis at the 95% 
and 99% confidence level are listed in Table [TJ all the cross section models implemented in the 



simulation exhibit equivalent behaviour regards the compatibility with the experimental data, 
with the exception of the Kahoul et al. model. The contingency table comparing the Kahoul 
et al. and ECPSSR-HS models confirms that the two models show a statistically significant 
difi^erence regards their accuracy (p-value of 0.001). 

The contingency tables associated with the other models show that they are statistically 
equivalent regards their accuracy. However, when only the lower energy range (below 5-7 MeV, 
depending on the atomic number) is considered, a statistically significant difference at the 95% 
confidence level (p-value of 0.034) is observed between the ECPSSR model and the ECPSSR-HS 
one; the latter is more accurate with respect to experimental data. 

From this analysis one can conclude that the implemented K shell ionization cross section 
models exhibit a satisfactory accuracy with respect to experimental measurements. 

L shell ionization cross sections 
The cross sections for L sub-shell ionization cross sections were compared to the experimental 
data collected in two complementary compilations [27j, [28]. The same method was applied as 
described for the validation of K shell cross sections. 

The ECPSSR model appears to provide a satisfactory representation of L shell ionisation 
cross sections with respect to experimental data, especially with its United Atom variant. 

The ECPSSR-UA exhibits the best overall accuracy among the various models; the Orlic et 
al. model exhibits the worst accuracy with respect to experimental data. This semi-empirical 
model is the only option implemented in Geant4 9.2 for the calculation of L shell ionization 
cross sections. 

The accuracy of the various cross section models was studied by means of contingency tables 
to evaluate their differences quantitatively. The categorical analysis was performed between 
the ECPSSR model with United Atom correction, i.e. the model showing the best accuracy 
according to the results of the test, and the other cross section models. The contingency 
tables were built based on the results of the test at the 95% confidence level, summing the 
"pass" and "fail" outcome over the three sub-shells. 

The Orlic et al. semi-empirical model is found to be significantly less accurate than the 
ECPSSR-UA model: the hypothesis of equivalence of their accuracy with respect to experimental 
data is rejected at 99% confidence level. The p-values concerning the comparison of the 
Miyagawa et al. empirical model are close to the critical region for 95% confidence, and slightly 
different for the three tests performed on the related contingency table. The Sow et al. empirical 
model and the ECPSSR model in its original formulation appears statistically equivalent in 
accuracy to the ECPSSR model with United Atom correction. 

As a result of this analysis, the ECPSSR model with United Atom approximation can be 
recommended for usage in Geant4-based simulation applications as the most accurate option 
for L shell ionization cross sections. The ECPSSR model in its original formulation can be 
considered a satisfactory alternative; the Sow et al. empirical model has satisfactory accuracy, 
but limited applicability regards the target elements and proton energies it can handle. 

4-3. Cross section models for high energy PIXE 

PIXE as a technique for elemental analysis is usually performed with proton beams of a few 
MeV. In the recent years, higher energy proton beams of a few tens MeV have been effectively 
exploited too. High energy protons are a source of PIXE in the space radiation environment. 

The interest in high energy PIXE has motivated recent theoretical investigations [TO] 
concerning cross section calculations at higher energies. Despite the emerging interest of high 
energy PIXE, only a limited set of experimental data is available above the energy range of 
conventional PIXE techniques. 



Table 2. P-values from the test concerning high energy experimental data 
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0.235 
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0.351 
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The accuracy of the implemented K shell cross section models was evaluated against two sets 
of measurements at higher energy |29l I30j . respectively at 66 and 68 MeV. The experimental 
measurement with uranium was not included in the comparison, since it appears affected by 
some experimental systematics. 

The -yd test was performed first separately on either experimental data set to evaluate the 
possible presence of any systematics in the two test cases, then on the combined data set. The 
p-values from the test against these experimental data are listed in Table [2l 

Over the limited data sample considered in this test, the ECPSSR model with the correction in 
[19Uhis model does not appear to provide better accuracy than the original ECPSSR formulation; 
nevertheless more high energy experimental data would be required to reach a firm conclusion. 
Also, this analysis should be verified over tabulation deriving from a published version of the 
ISICS code, when it becomes available. 

5. Application of the PIXE prototype software 

The prototype components for PIXE simulation described in the previous sections were applied 
to a study of the passive, graded Z shielding of the X-ray detectors of the eROSITA telescope 
|31j on the upcoming Russian Spectrum-X-Gamma space mission. 

The background spectra due to cosmic-ray protons were simulated for the three different 
eROSITA graded Z shield designs. A comparison of the results is depicted in Fig. [5l 

This application demonstrates that the developed software is capable of supporting concrete 
experimental studies. Nevertheless, the concerns outlined in the previous sections should be 
kept in mind: while the present PIXE simulation component can provide valuable information 
in terms of relative fiuorescence yields from inner shells, the intrinsic limitations of the mixed 
transport scheme in which ionization is modelled and the lack of cross section calculations for 
outer shells prevent an analysis of the simulation outcome in absolute terms. 

6. Conclusion and outlook 

This paper presents a brief overview of the status, open issues and recent developments of PIXE 
simulation with Geant4; a more extensive report of the underlying concepts, developments and 
results is available in [3]. 

The analysis of the related models currently available in Geant4 showed their limitations 
under various aspects: physics functionality, correctness of implementation, accuracy, consistent 
behavior, software design, performance and usability. 

The new developments represent a significant step forward regards PIXE simulation with 
Geant4. They extend the capabilities of the toolkit by enabling the generation of PIXE 
associated with K, L and M shells for protons and a particles; for this purpose a variety of cross 
section models are provided. The adopted data-driven strategy and the software design improve 
the computational performance over previous Geant4 models. The validity of the implemented 
models has been quantitatively estimated with respect to experimental data. The results provide 
objective guidance for the optimal selection of simulation models in user applications. 




(a) Cu Shield 




(b) Cu-Al Shield 




(c) CU-AI-B4C Shield 

Figure 5. A comparison of the fluorescence background due to ionization by cosmic-ray protons 
in an L2 orbit for three different graded Z shield designs for the eROSITA X-ray detectors. 



An extensive ionisation cross section data library has been created as a by-product of the 
development process: it can be of interest to the experimental community for a variety of 
applications, not necessarily limited to PIXE simulation with Geant4. 

Some issues identified in the course of the development process are still open: they concern 
the consistency of PIXE simulation in a mixed condensed-discrete particle transport scheme. 
A possible approach to address them, while preserving the current continuous-discrete design 
scheme of hadron ionisation in Geant4, would involve the calculation of ionization cross sections 
for outer shells than M, This calculation is feasible, exploiting known theoretical methods, yet 
it would require a significant investment of resources. In parallel, a project [32] is in progress 
to address design issues concerning co-working condensed and discrete transport methods in a 



general purpose simulation system. 

Despite the known limitations related to mixed transport schemes, the software developments 
described in this paper provide sufficient functionality for realistic experimental investigations. 
The plans concerning the public availability of these new developments are discussed in |3j. 
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